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Abstract. - The solid-fluid phase transition of a granular material shaken horizontally is 
investigated numerically. We find that it is a second-order phase transition and propose two 
order parameters, namely the averaged kinetic energy and the averaged granular temperature, 
to determine the fluidization point precisely. It scales with the acceleration of the external 
vibration. Using this fluidization point as critical point, we discuss the scaling of the kinetic 
energy and show that the kinetic energy and the granular temperature show two different 
universal critical point exponents for a wide range of excitation amplitudes. 



Granular materials show a fascinating behaviour under a variety of experimental conditions O . 
Most striking are the different regimes found under vertical shaking: for low excitations sound 
waves are found, for medium excitations heap formations and convection rolls are visible 
and for high excitations period doubling can occur which can lead to standing and travelling 
waves |Q. This behaviour is connected to two visible phase transitions where the strength 
of the excitation is measured in the dimensionless variable T := Aoiv/g, with Aq and ui 
standing for the amplitude and the angular frequency of the external shaking and g denoting 
the gravitational constant. The first transition occurs around T — 1 ||, and the second was 
found to occur around T = 1.8 in a one-layer system In a quasi two-dimensional box, the 
transition points were significantly higher In the intermediate regime, the bed expansion 
scales as T M , whereas the supercritical scaling is correlated with Aq ui [|[ ^, ^, ^] . In this 
letter, we will focus on a slightly different setup and study the first phase transition from a 
solid-like to a fluid-like regime of a granular material undergoing horizontal vibrations. Our 
aim will be to quantify this transition point and to show the correct scaling with Aq and u>. 
Furthermore, we demonstrate that the averaged kinetic energy and the granular temperature 
show beautiful scaling over more than one decade with different critical exponents. 

The system in mind was first described by Strafiburger et al. || and an enhanced setup 
was used in J^] for quantitative measurements which we will sketch briefly. A few layers of 
commercially available glass ballotonies were filled into a 100 mm long and 0.6 mm wide 
channel and exposed to horizontal vibrations of the form A(t) = Aq sin(wt) where u = 2nf. In 
order to get better defined surface properties for particle- wall contacts a monolayer of smaller 
spheres was glued onto the channel bottom. Since the channel width is equal to the maximum 
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particle diameter, one gets a quasi two-dimensional system which can easily be observed from 
the side. For low excitations, cither low frequencies / or small amplitudes Aq, no relative 
motion of particles with respect to the container is visible and the spheres form a nearly 
perfect triangular lattice, sketched in fig. nk. For higher excitations, relative motion sets in 
and the lattice breaks up, sketched in fig. Tp. Using this system, the phase diagram for the 
solid-fluid transition was investigated and it was found that the transition point scales with T 
using experiments and numerical simulations J9| . A summary of the phase diagram is shown 
in fig. demonstrating that the numerical model accurately describes the experiment and the 
slight shift could be attributed to the higher wall friction present in the experiment. For the 
experimental values, the onset of grain motion is denoted by f and the extrapolated critical 
frequency by f c . 
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Fig. 1. - Sketch of the physical system undergoing horizontal vibrations for (a) low and (b) higher 
excitations. 



Fig. 2. - Phase diagram for the solid-fluid transition given in ref. M: experimental points are * (fo) 
and • (fc); data from numerical simulations are o with solid line as best fit. Dashed line shows 
theoretical curve for equal size particles. 



We will study the above mentioned system numerically using spherical particles which are 
bound to move in two dimensions and only interact via contact forces. Periodic boundary 
conditions are used in the direction of the shaking (horizontal direction) and 342 particles are 
placed to form two layers initially. A predictor-corrector time integration scheme commonly 
used in molecular dynamics type simulations is used. The forces acting on particle i, having 
a radius of r» = d%/2, during a collision with particle j are in the normal direction (n) 

F'ij = -kn(n + rj - njfi) - 7„%n (1) 

and in the shear direction (s) 

F? d = sign(uy-s) • min(7 s *%s, fJ,\F^\) . (2) 
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Here stands for the vector joining both centers of mass and Vij denotes the relative motion of 
the two particles. The particles are allowed to rotate since this is well observed experimentally. 
The parameter k n is related to the Young modulus and is chosen as 10 6 N/m which is high 
enough to avoid unphysical results due to the detachment effect JTo[ ] . The parameter j n is 
related to the experimentally more accessible restitution coefficient e n and a value of e n = 0.75 
was used in all simulations in accordance with the experiment ||. The relevant parameters 
entering the shear force are given by the detailed collision experiments conducted by Foerster 
et al. (II]]. Our parameter 7 S was set to such a high value that the shear force was mostly 
dominated by the Coulomb threshold condition [i\F™-\ where a value of /i = 0.1 was used for 
particle-particle and [i w — 0.13 for particle- wall contacts. Our force law in the shear direction, 
eq. (P), is a slight oversimplification with respect to the experimental results in the low shear 
11 since it does not contain the reversion of rotation. We checked the validity 
of this approximation by using a shear force that exactly matches the data points given in 
ref. JllJ for glass spheres and even though the transition point shifted slightly neither the 
scaling of the transition points with T nor the critical point exponents changed their values 
within the numerical accuracy. A more detailed discussion of the different force laws was given 
by Schafer et al. E3. A system size of 10.08 cm was used and the average particle diameter 
was (d) = 0.56 mm with a uniform size variation of ± 0.04 mm. The container bottom consists 
of smaller particles with an average size of (d) / 4 denoted by black particles in fig. [j]. 

In order to detect the relative motion between particles, we look at the average kinetic 
energy of all particles, defined as ej tm (r) := ^ Y2i=i m i v i( T )> as function of time measured 
in the moving reference frame of the channel. This signal is averaged over an integer number 
of external vibration cycles, typically 5-10, to give E^ n and we propose that it can be used 
as an order parameter for the solid-fluid transition in granular assemblies. In fig. |^a, we 
show how E^ n depends on the shaking frequency / measured in Hz, for two different shaking 
amplitudes Aq (• - 7 cm, o - 3 cm). The errorbars are less than the symbol size. Below a 
well defined threshold / c , £^i n is approximately zero which indicates that the particles are all 
at rest and can be viewed as a solid. Above the threshold, a monotonic, drastic increase is 
visible. Particles are moving around relative to each other and the average density decreases. 
The state of the system in this regime can be best described by a fluid. 

A more accurate method to determine f c is given by subtracting the bulk motion from 
the kinetic energy in the following fashion tgran (V) := ^iv Eili m i{ v i{ T ) ~ ( v )) 2 which is 
commonly referred to as the granular temperature of the system. Here (v) stands for the 
average particle velocity at time r. This quantity is also averaged over an integer number of 
full cycles of the external vibration to give Tg ran . Below the transition point, this quantity is 
zero and it is non-zero above it. Since Tgran scales linearly with / for a wide frequency range 
it can easily be extrapolated to the limit Tgran — > to give a very accurate value of f c . This 
procedure is demonstrated in fig. ||b for the two shaking amplitudes Aq = 7 cm and Aq = 3 
cm and gives values of / c = 0.91 and 1.4 Hz, respectively. When kinetic theories were applied 
to vibrated beds the granular temperature was found to scale as T ~ (A uj) 2 || 0] which is 
in contrast to our finding and indicates that the current theories do not adequately describe 
the solid-fluid transition in granular materials. 

When the critical frequencies for an amplitude range of 0.5 to 9 cm are plotted on a double 
logarithmic scale, one finds a scaling of Aq/ 2 = const, which is shown in fig. ^. Such a 
scaling relation was also found for the transition point for surface instabilities and the onset 
of convection rolls in a sand pile under vertical vibrations || . It can be understood by looking 
at the force balance along the shear direction of one particle on a row of particles. This gives 
as necessary condition for particle motion Aq uj 2 / g > tan© where O measures the angle from 
the vertical of the center of mass of the top particle to the next particle center below it. For 
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Fig. 3. - Averaged kinetic energy (a) and granular temperature (b) as function of frequency / for two 
different shaking amplitudes A) = 7 cm and Ao = 3 cm. 



equal size spheres this gives O e = 30° whereas in our case 6 = 11.5°. 

Knowing that the transition point scales as Aof 2 = const., we replot in fig. [| the av- 
eraged kinetic energy -E^-in divided by Aq as function of the dimensionless quantity T' := 
A-of 2 /(g tan 9). All data points collapse onto one straight line, giving a transition point of 
T' c k 1.07. For the vertical shaking experiments, a transition point of r c w 1.2 was found M. 

The transition point does not depend on the choice of the restitution coefficient in the normal 
direction, even though it changes the kinetic energy and has a dramatic effect on the absolute 
value of the granular temperature. On the other hand, changing the Coulomb threshold /j. 
shifts the critical point when a high enough value of 7 S is used. For /j, — 0, the transition 
frequency is lowered by roughly 0.5 Hz; whereas increasing it to [i = 1 increases f c by up to 
1 Hz as well. Changing the system size to 7.56 or 15.12 cm did neither change the transition 
point nor the absolute values of the averaged kinetic energy and granular temperature. 
We will now turn to the main issue of this letter, namely the calculation of critical-point 
exponents for the solid-fluid phase transition in granular material. The sharp transition at 
f c from a zero value kinetic energy (or granular temperature) to a monotonically increasing 
value with increasing frequency suggests that the solid-fluid phase transition in a granular 
material undergoing vertical vibrations can be described by a second-order phase transition. 
This is stressed by the fact that no hysteresis can be found in the numerical simulations when 
increasing or decreasing the driving frequency, /. It was shown above that the kinetic energy 
or the granular temperature can both equally well be used as order parameter to characterize 
this transition. We use 

f-fc 



' fc 

as dimensionless variable to describe the behavior near the critical point f c 
exponent A is then given by 



131. The critical 



A:=lim^ 

e-^o In e 
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Fig. 4. - Critical frequency as function of the external shaking amplitude Ao 

Fig. 5. - Rescaled averaged kinetic energy as function of dimensionless quantity F' := Aof 2 /(g tan G). 

where h{e) stands for either the kinetic energy or the granular temperature Tg ra n- 

Assuming the existence of this limit, we can approximate h(e) in the vicinity above the critical 
point by 

h(e) ~ 6 A . 

We applied this technique to 13 different shaking amplitudes Aq in the range 0.5 to 9.0 cm 
and found an universal exponent for the scaling of the kinetic energy of 1.25 ± 0.15 and of the 
granular temperature of 1.0 ± 0.15. The error bars stem from the possible error in the critical 
frequency. But no systematic trend of the exponent for decreasing or increasing values of Aq 
can be observed. In fig. ||, we show this scaling for shaking amplitudes of Ao = 7, 3 and 1 cm 
using either the kinetic energy and an exponent of 1.25 (a) or the granular temperature and 
an exponent of 1.0 (b) for all three curves. For each curve, we used values of e ranging over 
at least one order of magnitude. For Aq — 3 cm, we extended our data points to very high 
excitation amplitudes and find a different exponent for values of e > 3 which might indicate 
that the transition from the fluid to a more gas-like state takes place in this regime || . For the 
lowest shaking amplitude, the left most point shows a slight deviation from the fit which can 
be attributed to a less accurate value for the transition frequency for this excitation amplitude 
and finite size effects. 

In this letter, we investigated numerically the solid-fluid phase transition of a granular 
material exposed to horizontal shaking. By looking at the averaged kinetic energy or the 
granular temperature of all particles, we found that these quantities are only non-zero above a 
well defined transition frequency which we called critical frequency and showed no hysteresis. 
This suggests that this phase transition is of second-order and we presented the correct scaling 
of the transition frequency with respect to the shaking amplitude. In the vicinity of the 
critical point, we calculated the critical exponents associated with the kinetic energy and the 
granular temperature and found an universal exponent of 1.25 for the first and of 1.0 for the 
latter quantity. The same technique can probably be applied to the solid-fluid phase transition 
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under vertical vibration and it would be interesting to compare the exponents with our results. 

*** 

We would like to thank I. Rehberg, G. StraBburgcr, S. Grossmann, A. Esser and K. Drese for 
enlightening discussions. This work was supported by the Deutsche Forschungsgemeinschaft 
through Ri 826/1-1. 

REFERENCES 

[1] H.M. Jaeger, S.R. Nagel and R.P. Behringer, Physics Today, 4 (1996) 32 
[2] S. Warr, J.M. Huntley and G.T.H. Jacques, Phys. Rev. E, 52 (1995) 5583 and references 
cited therein 

[3] P. Evesque and J. Rajchenbach,, Phys. Rev. Lett., 62 (1989) 44; E. Clement, J. Duran 

and J. Rajchenbach,, Phys. Rev. Lett, 69 (1992) 1189 
[4] C.E. Brennen, S. Ghosh and C. Wassgren, in Powders & Grains 93„ edited by C. Thornton, 

(Balkema) 1993 

[5] A. Goldshtein, M. Shapiro, L. Moldavsky and M. Fichman, J. Fluid Mech., 287 (1995) 
349 

[6] S. Luding, H.J. Herrmann and A. Blumen, Phys. Rev. E, 50 (1994) 3100 
[7] J. Lee, Physica A, 219 (1995) 305 

[8] G. Strassburger, A. Betat, M.A. Scherer and I. Rehberg, in Traffic and Granular Flow,, 
edited by D.E. Wolf, M. Schreckenberg and A. Bachem, (World Scientific) 1996 

[9] G.H. Ristow, G. Strassburger and I. Rehberg, Phys. Rev. Lett., 79 (1997) 833 
[10] S. Luding, E. Clement, A. Blumen, J. Rajchenbach and J. Duran, Phys. Rev. E, 50 (1994) 
4113 

[11] S.F. Foerster, M.Y. LOUGE, H. Chang and K. Allia, Phys. Fluids, 6 (1994) 1108 
[12] J. Schafer, S. Dippel and D.E. Wolf, J. Phys. I France, 6 (1996) 5 

[13] H.E. Stanley, Introduction to Phase Transitions and Critical Phenomena, (Clarendon Press) 
1971 



